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摘要 : 以 印度 块 菌 (Tuber indicum). 丰产 区 云南 省 为 例 ， 利 用 地 理 信息 系统 和 物种 分 布 模型 ， 并 通过 影响 
印度 块 菌 分 布 的 主要 环境 因子 和 和 气候 因子 ， 共 同 模拟 现在 及 未 来 生境 分 布 模式 。 结 果 表 明 结 合 环境 因子 和 
气候 因子 ， 在 建 模 过 程 中 能 提高 模型 预测 准确 度 ， 在 几 类 物种 分 布 模型 中 ，MAXENT 模型 具有 最 优 的 拟 
合 效 果 。 在 大 的 空间 尺度 上 ， 年 降水 ， 最 湿 季 度 降水 ， 最 冷 月 份 最 低温 、 地 貌 类 型 及 土壤 类 型 对 印度 块 菌 
的 生境 分 布 影响 最 大 。 此 外 ， 在 模型 模拟 的 A2 和 B2 未 来 气候 变化 情景 下 ， 印 度 块 菌 未 来 的 新 增生 境 均 
呈 北 上 趋势 ， 且 B2 情景 下 生境 的 适 生 程度 低 于 A2 情景 。 
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Abstract: In order to facilitate conservation of forest products, a case study focused on monitoring the distribution 
trend of an economic mushroom, Tuber indicum, was conducted in Yunnnan province. Geographic Information Sys- 
tem and species distribution models were employed to establish both current potential habitat and a future projection 
of T. indicum habitat, through main environmental and climate variables. Results confirmed that combining environ- 
mental and climatic factors in the model building process could greatly increase model accuracy. Among tested spe- 
cies distribution models, MAXENT showed the highest performance. At a large geo-spatial scale, annual precipitati- 
on during the wettest quarter, minimum temperature during the coldest month, geomorphology and soil types played 
the most significant roles in determining T. indicum distribution. Our future projection showed T. indicum shifting 
northward under both modeled climate change scenarios ( A2 and B2) , with the B2 scenario having a lower degree of 
habitat suitability. 
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印度 块 菌 (Tuber indicum) 是 一 种 重要 的 经 
济 食用 菌 ， 通 常 以 外 生 菌 根 形式 与 寄主 植物 形成 
共生 关系 ， 常 见于 海拔 2 000 ~3 000 m 的 针 叶 
林 、 阅 叶 林 或 混交 林 中 〈 陈 应 龙 和 号 明 钦 ，2000; 
Zhang 等 , 2005) 。 子 实体 生长 于 地 下 10 ~ 15 cm， 
成 熟 时 散发 的 芳香 类 物质 能 强烈 吸引 咕 齿 类 哺乳 
动物 ， 被 认为 具有 动物 性 激素 作用 ( 陈 惠 群 等 ， 
1998a; Gao 等 , 2004) 。 在 欧洲 ， 训 练 块 菌 犬 是 
寻找 此 类 地 下 真菌 的 主要 方式 (Claus 等 ， 
1981 ) 。 作 为 外 生 菌 根 真 菌 ， 印 度 块 菌 通过 产生 
菌 根 结构 对 共生 树种 产生 影响 ， 从 而 参与 到 整个 
生态 系统 循环 过 程 〈( 柳 成 益 等 ,2008; Geng 等 ， 
2009) 。 此 外 ， 印 度 块 菌 具有 良好 的 营养 保健 功 
能 〈 陶 恺 和 刘波 ,1990) ， 能 有 效 的 增强 机 体 免 
疫 能 力 并 预防 癌症 (HAWE, 1994), —E A 
来 受到 人 们 喜爱 。 

印度 块 菌 在 我 国 大 部 分 地 区 广泛 分 布 ， 据 报 
道 ， 云 南 和 四 川 两 省 是 其 主要 产 区 ( 任 德军 等 ， 
2005; García-Montero 等 ，2010)。 印 度 块 菌 是 当 
地 重要 的 林 下 产品 之 一 ， 对 经 济 发 展 具 有 推动 作 
用 。 然 而 ， 目 前 广泛 采用 的 原始 粗放 泛 猎 式 采集 
方式 已 造成 了 原生 生境 的 极 大 破坏 ， 印 度 块 菌 已 
归 为 高 等 真菌 类 濒危 类 群 ( 谭 著 明和 付 绍 春 ， 
2002; 刘 培 贵 等 , 2003; 王 云 和 刘 培 贵 ，2009 ) 。 
为 了 应 对 日 益 增 长 的 市 场 需求 ， 并 对 森林 自然 资 
源 进行 合理 保护 ， 中 国 自 20 tita 80 年 代 后 期 便 
开始 探索 人 工 栽培 块 菌 技 术 ， 并 于 1989 年 建立 
了 最 早 的 块 菌 种 植 园 ( Hu 等 , 2005) ， 随 后 ， 在 
贵州 、 云 南 、 四 川 、 湖 南 和 广东 省 也 开展 了 相关 
人 研究 (García-Montero 等 , 2010 ) 。 人 工 栽 培 块 菌 
的 成 功 ， 极 大 缓解 了 市 场 需求 压力 ， 也 扩大 了 印 
度 块 菌 的 原生 生境 范围 〈 陈 娟 等 , 2009 ) 。 

国内 外 学 者 对 块 菌 的 研究 多 基于 生物 因素 ， 
如 寄主 植物 、 小 型 哺乳 动物 、 土 壤 动 物 、 土 壤 微 
生物 、 人 类 因素 等 (Lehmkuhl 55, 2004; 苏 开 
美 , 2007); 非 生物 因素 ， 如 土壤 理化 性 质 、 降 
水 、 温 度 、 地 势 等 对 块 菌 菌 丝 生 长 、 子 实体 大 小 、 
成 熟 周期 、 产 量 情况 等 影响 (Tang 等 , 2008) 却 
很 少 和 生物 因素 结合 起 来 ， 共 同 预 测 块 菌 的 适宜 
生境 。 现 有 的 生境 适宜 性 评估 是 基于 传统 的 生态 
学 方法 之 上 的 ， 如 唐 平 等 (2005) 采用 Vander- 
loeg 和 Scavia 选择 系数 Wi 和 选择 指数 Ei 作为 评 






































价 标准 ， 对 攀枝花 地 区 的 块 菌 适 宜生 境 进 行 了 探 
究 。 陈 慧 群 等 (1998b) 通过 对 块 菌 原 产 地 的 气 
候 分 析 ， 提 出 了 地 温和 降雨 量 是 对 块 菌 生长 影响 
较 大 的 气象 因素 ， 然 而 ， 这 些 研 究 主要 基于 局 部 
地 区 监测 ,缺乏 时 间 和 空间 尺度 上 的 动态 模拟 ， 
难以 在 其 它 地 区 通用 。 

全 球 气 候 变 暖 已 是 不 可 忽视 的 世界 性 话题 ， 
中 国 近 100 年 来 平均 地 表 温 度 上 升 了 0.5 ~0.8Y， 
增 温 速率 较 全 球 其它 地 区 略 强 ( 丁 一 汇 等 ， 
2009) 。 已 有 研究 表明 气候 变化 已 经 对 物种 分 布 
造成 影响 ( Parmesan 和 Yehe, 2003; Root 等 ， 
2003), ， 未 来 可 能 发 挥 更 大 作用 ( Peterson 等 ， 
2001; Beaumont 和 Hughes, 2002; Araájo 等 ， 
2004; Thuiller, 2004; Zonneveld £, 2009; Saran 
等 , 2010; Luedeling 等 , 2011) 。 不 同 物种 对 气候 
变化 的 响应 存在 差异 ， 根 据 观测 资料 显示 ， 全 球 
尺度 上 的 变 暧 可 能 会 导致 大 型 真菌 子 实体 生长 
季节 提前 及 生长 期 延长 等 现象 (Kauserud 等 ， 
2008; Kauserud 等 ,2010 ) 。 近 期 研究 还 发 现 除 
了 对 子 实 体 生 长 有 影响 ， 未 来 气候 变化 还 可 能 
SUM EBZXERS (Büntgen SF, 2011) 。 本 文 基于 
地 理 信息 系统 和 不 同 的 物种 分 布 模型 ， 首 先 通过 
影响 印度 块 菌 分 布 的 主要 环境 因子 和 气象 因子 ， 
对 云南 省 印度 块 菌 生境 分 布 现状 进行 模拟 , 采用 
曲线 下 面积 (Area Under Curve, 缩写 为 AUC) 
筛选 最 优 模型 及 参数 ， 并 用 于 建立 未 来 不 同 气候 
情景 下 的 印度 块 菌 分 布 模型 。 考 虑 到 印度 块 菌 与 
共生 树种 的 紧密 联系 ， 建 立 分 布 现状 模型 时 加 入 
了 植被 因子 ， 提 高 了 模型 准确 性 。 而 植被 本 身 对 
气候 变化 也 有 响应 (Loehle 和 Blanc, 1996) ， 故 
本 研究 采用 了 保守 的 建 模 方式 来 分 析 未 来 气候 变 
化 情景 下 的 印度 块 郑 空间 分 布 情况 ， 噜 除了 植被 
因子 来 建立 模型 ， 在 一 定 程度 上 减少 了 模型 的 不 
确定 性 。 本 研究 旨 在 对 将 来 块 菌 种 植 园 选 址 提供 
参考 依据 ， 同 时 也 为 其 它 森 林 资 源 的 合理 管 护 黄 
定 基础 。 






































1 材料 和 方法 
1.1 研究 区 域 概况 

云南 省 位 于 ES97931'39"-106?11'47", N 21?08'32"— 
29°15'8" 之 间 ， 土 地 总 面积 为 39.4 万 km ， 占 全 国 陆地 
总 面积 的 4.1% 。 省 内 8096 为 山地 类 型 ， 地 势 北 高 南 低 ， 
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垂直 高 度 为 84 ~6 044 m (图 1) 。 为 典型 低 纬度 高 海拔 
地 区 ， 具 有 四 季 温 差 小 , 干 湿 季 节 分 明 的 特点 。 由 于 同 
时 受到 西南 面 印度 洋 、 东 南面 太平 洋 季 风 及 西北 面 青藏 
高 原 气候 的 影响 ， 云 南 形 成 了 独特 的 水 平地 带 性 和 垂直 
地 带 性 。 云 南 省 是 全 国 植物 区 系 重要 组 成 部 分 ， 主 要 
植被 类 型 包括 热带 雨林 、 热 带 季 雨林 、 和 常 绿 阔 叶 林 、 
落叶 闪 叶 林 、 硬 叶 常 绿 阔 叶 林 、 针 叶 林 等 ， 为 野生 动物 
提供 了 丰富 的 栖息 地 (Parmesan 和 Yehe, 2003; Brown 
和 Davis, 2006) 。 全 省 野生 食用 菌 约 750 多 种 ， 占 全 国 
时 生食 用 菌 种 类 88.24% ， 其 中 具有 野生 贸易 真菌 64 属 
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图 1 云南 省 海拔 梯度 图 


Fig.1 The elevation gradients in Yunnan Province 


207 种 ( 王 向 华 和 刘 培 贵 , 2002; 丁 一 汇 等 , 2009)， 故 
素来 享有 “植物 王国 "”、“ 动 物 王国 ”和 “食用 菌 王国 ” 
的 美誉 ， 是 全 球 生物 多 样 性 热点 地 区 之 一 ( 苏 开 美 ， 
2005), 

1.0 研究 方法 

1.2.1 印度 块 菌 分 布 数据 收集 通过 查阅 昆明 植物 研 
究 所 隐 花 植物 标本 馆 (HKAS) 印度 块 菌 历 史 采 集 记 录 ， 
参考 相关 文献 和 报告 ， 并 结合 野外 实际 调查 共同 收集 分 
布 数据 。 采 用 ArcGIS 的 空间 统计 分 析 工 具 (Spatial Sta- 
tistics Tools), Alf 1 km 内 重复 记录 、 空 间 自 相关 及 人 
工种 植 块 菌 地 点 后 ， 最 终 保 留 46 个 分 布点 。 从 国家 基础 
地 理 信 息 系 统 (http://nfgis. nsdi. gov. cn/) 下 载 1: 400 
万 全 国 矢量 地 图 作为 底 图 ， 并 将 获取 到 的 坐标 信息 保存 
于 Excel 2007 电子 表格 。 
1.2.2 气候 数据 及 环境 因子 收集 从 “WorldClim” 数 
据 库 (http://www. worldclim. org/bioclim. htm) FRAIR 
数据 ， 分 辩 率 为 1 km。 该 数据 库 覆 盖 了 全 球 表面 (土地 
面积 不 包括 南极 洲 ) 大 部 分 地 区 ， 是 Robert 等 通过 收集 
全 球 各 地 1950-2000 年 气象 台 站 月 降水 量 、 降 雨量 和 温 
度数 据 集 ， 采 用 ANUSPLIN 数据 包 中 的 thin-plate smoot- 
hing spline 算法 推导 后 建立 。 以 往 对 块 菌 的 生态 地 理 研 
究 表明 ， 除 气候 外 〈 陈 惠 群 等 ,1998b) ， 植 被 、 土 壤 、 
地 形 因素 等 也 对 其 分 布 起 重要 作用 ( 陈 惠 群 等 ，1999; 
Trappe 等 , 2009) 。 如 在 中 国 西南 地 区 ， 印 度 块 菌 多 发 生 
于 1000 ~3000m 山区 ， 坡 度 通常 小 于 35。 的 阴 坡 或 半 阴 
坡 生 境 ， 其 主要 寄主 植物 包括 云南 松 (Pinus yunnanen- 
sis) 、 华 山 松 (P. armandi) 、 栎 树 (Quercus spp. ) , T&K 
(Alnus nepalensis) 、 云 南 油 杉  ( Keteleeria evelyniana) 、 水 
ZEX (Viburnum cylindricum). 等 (WW JE, 2000; 苏 开 
X, 2005; Long 和 Li, 2009) ， 并 主要 生长 于 呈 弱 碱 性 或 
中 性 的 紫色 土 和 石灰 岩 土壤 中 (张大 成 和 王 云 , 1990), 
最 终 ， 根 据 数据 的 可 获取 性 ， 选 择 了 6 个 环境 因子 图 层 
用 于 建 模 ( 表 1)。 





















































































































































表 1 环境 因子 属性 列表 


Table 1 List for environmental factors attribute 
































» i 尽 / 分 辩 率 39r HEA H 
环境 因子 Factor 来 源 Source 比例 A eR 数据 类 型 

Scale/Resolution Data type 
植被 类 型 云南 植被 编写 组 , 1987 M o0 类 别 变量 
Vegetation types Yunnan vegetation compiling group, 1987 ; Categorical 
土壤 类 型 云南 省 国土 资源 地 图 集 编 委员 会 , 1990 eooo O00 类 别 变量 
Soil types Yunnan province national land and resource atlas compiling group, 1990 i Categorical 
HKA 云南 省 国土 资源 地 图 集 编 委员 会 , 1990 boo 类 别 变量 
Geomorphology Yunnan province national land and resource atlas compiling group, 1990 : Categorical 
坡度 、 坡 向 和 海拔 SRTM 数字 高 程 模型 90 m 连续 变量 
Slope, aspect and altitude Continuous 


SRTM Digital Elevation Model 
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1.2.3. 未 来 气候 研究 ”将 1950-2000 年 气候 设 定 为 基 
准 线 (Baseline)， 从 气候 包 络 (Climate Envelope) 中 可 
直接 获取 19 个 生物 气候 变量 ， 而 未 来 的 3 个 时 段 ， 即 
Baseline-2020, 2020—2050 , 2050-2080 年 的 生物 气候 变 
量 需 通过 DIVA-GIS 软件 从 最 高 温 、 最 低温 、 降 雨 数据 
和 海拔 数据 共同 推导 得 出 。 本 研究 采用 CCCMA ( Cana- 
dian Centre for Climate Modelling and Analysis ) 全 球 气候 
模型 模拟 的 A2 和 B2 未 来 气候 情景 。 根 据 政府 间 气 候 变 
化 委员 会 (IPCC) 排放 情景 特别 报告 (SRES) ， 在 A2 
气候 情景 下 全 球 各 部 分 变化 存在 高 度 异 质 性 。 主 要 特点 
是 人 口 持续 增长 ， 经 济 发 展 较 其 它 情 景 模式 缓慢 ， 常 依 
赖 于 局 部 或 区 域 性 资源 分 布 ， 发 展 速度 不 均一 且 资 源 的 
利用 较为 分 散 ， 注 重 区 域 间 合 作 。B2 气候 情景 描述 了 区 
域 性 经 济 、 社 会 和 环境 可 持续 利用 发 展 模式 ， 人 口 增长 
较 A2 情景 缓慢 ， 经 济 发 展 速度 居中 ， 具 有 逐步 改善 的 
多 元 化 新 技术 ， 更 注重 环境 保护 和 社会 公平 Naki- 
cenovic 等 , 2000) 。 

1.2.4 建 模 方 法 “在 各 种 物种 分 布 模型 中 ， 最 大 精 模 
型 (Maximum Entropy Method, 缩写 为 MAXENT) 是 一 种 
较 新 的 算法 ， 具 有 其 特有 的 优越 性 ， 能 够 在 最 小 样本 的 
情况 下 得 到 较 好 的 模拟 结果 (Graham 和 Hijmans, 2006; 
Hernandez 等 , 2006), 。MAXENT 模型 已 经 成 功 运用 于 模 
拟 现 在 和 将 来 植物 的 分 布 ， 如 对 Banksia priontes 的 生境 
模拟 (Elith 等 , 2011) 。 预 测 的 潜在 适宜 地 点 需 能 最 大 
程度 满足 研究 物种 对 自然 环境 的 需求 条 件 (Phillips 和 
Dudik, 2008) 。 同 时 ， 本 研究 中 采用 的 BIOCLIM 和 DO- 
MAIN 生态 位 模型 是 内 髋 于 DIVA-GIS 中 的 气候 模型 , 与 
MAXENT 模型 相 类 似 ， 能 够 利用 不 完整 的 信息 来 预测 物 
种 潜在 栖息 地 (Carpenter 55, 1993; Beaumont 等 ,2005 ) 。 
只 需要 物种 存在 资料 和 气候 包 络 ， 便 能 对 潜在 和 未 来 的 
生境 区 域 进 行 预测 。 背 景 像 元 被 用 来 代替 物种 在 真实 环 
境 中 的 不 存在 数据 ， 能 够 避免 通常 情况 下 物种 存在 或 不 
存在 (Presence/Absence) 模型 中 容易 导致 的 假 阳性 错误 。 
1.2.5 数据 预 处 理 ” 需 将 各 个 环境 因子 图 层 转 换 为 . grd 
格式 (BIOCLIM 和 DOMAIN) 或 . asc 格式 ( MAXENT)， 
在 建 模 前 ， 利 用 ArcGIS 的 空间 主 成 分 分 析 功 能 ( Princi- 
ple Component Analysis) 去 除 相 关 性 大 于 0.9 的 图 层 ， 以 
避免 由 于 图 层 间 的 空间 高 度 相 关 性 带 来 的 模型 噪音 。 将 
样 点 数据 随机 分 为 两 部 分 ，50% 为 训练 数据 集 ，50% 为 
检测 数据 集 。 

1.2.6 模型 评估 MAXENT，BIOCLIM 和 DOMAIN 将 
同时 用 于 建立 1950-2000 年 印度 块 菌 潜在 分 布 图 ， 相 对 
于 以 往 只 考虑 气候 因子 的 模型 ， 本 研究 还 考虑 了 物种 之 
间 相 互 关 系 对 模型 的 影响 ，( Aratijo 和 Luoto, 2007) ， 
此 ， 共 有 三 类 模型 组 成 : (1) 气候 模型 ， 即 只 采用 独立 
的 气候 因子 用 于 建 模 ; (2) 气候 物种 模型 ， 即 除了 气候 




















































































































































































































因子 ， 还 包括 了 表 1 列举 的 所 有 环境 因子 ; (3) 气候 地 
形 模型 ， 即 将 气候 物种 模型 中 的 植被 类 型 剔除 后 的 情 
况 ， 只 采用 相对 稳定 的 地 形 、 土 壤 等 因子 用 于 建 模 ， 是 
一 种 保守 的 预测 方式 。 保 证 模型 输入 参数 一 致 ， 每 种 模 
型 运行 500 次 ， 随 机 选取 10 组 运行 结果 ， 计 算 曲 线 下 面 
积 的 平均 值 和 标准 差 ， 用 于 三 种 模型 间 预 测 准确 度 比 
较 。 选 择 准 确 度 最 高 的 模型 ， 分 析 各 环境 因子 和 气候 因 
子 对 建 模 的 重要 性 程度 ， 对 模型 贡献 的 累计 频率 达到 
90% 以 上 的 因子 将 被 视 为 主导 因子 。 

1.2.7 未 来 气候 变化 情景 建 模 ”理论 上 而 言 ， 要 预测 
印度 块 菌 分 布 对 未 来 气候 变化 的 响应 ， 必 须 掌 握 至 少 两 
方面 的 信息 : (1) 影响 印度 块 菌 分 布 的 环境 因子 在 未 来 
的 变化 情况 ; (2) 印度 块 菌 与 未 来 环境 因子 之 间 的 相互 
关系 。 而 事实 上 ， 目 前 上 述 信息 是 未 知 的 ， 只 能 根据 人 
们 对 过 去 和 现在 的 情况 的 认 知 来 进行 预测 ， 并 基于 该 预 
测 来 进一步 推测 物种 分 布 的 可 能 变化 。 因 此 ， 预 测 是 在 
特定 条 件 下 的 可 能 发 生 情况 。 本 研究 探讨 的 是 未 来 气候 
变化 情景 下 印度 块 菌 的 分 布 变化 ， 为 了 降低 模型 的 不 确 
定性 ， 假 定 其 它 条 件 保持 不 变 。 本 研究 主要 有 三 个 基本 
假设 : (1) 影响 印度 块 菌 分 布 的 环境 因子 中 ， 假 定 非 生 
物 环 境 因 子 ， 如 土壤 、 地 形 、 地 貌 等 在 未 来 保持 相对 恒 
定 ; (2) 印度 块 菌 与 现 有 寄主 树种 之 间 的 共生 关系 保持 
不 变 ; (3) 印度 块 菌 当代 和 未 来 的 生态 位 保持 不 变 。 为 
了 便于 比较 未 来 印度 块 菌 的 生境 范围 变化 情况 ， 首 先 将 
基准 线 图 和 不 同时 段 的 预测 图 以 二 值 图 层 (0 和 1) JÉ 
式 输出 ,0 代表 非 生境 区 域 ，!1 代表 生境 区 域 , 通过 DI- 
VA-GIS 将 未 来 生境 图 与 基准 线 图 重 且 (overlay)，, 采用 
图 层 相 减 (subtract) 功能 获取 三 个 时 段 ， 即 1950 - 
2020, 2020-2050 和 2050-2080 的 生境 迁移 图 。 将 出 现 
的 四 类 情况 分 别 是 : 未 来 减少 的 生境 (01); 非 适 生 区 
(00); 未 来 无 变化 生境 (11) 及 未 来 新 增生 境 (10). 
为 了 进一步 揭示 生境 变化 过 程 中 适 生 程度 的 变化 ， 生 成 
未 来 不 同 气候 情景 下 印度 块 菌 分 布 的 概率 图 ， 根 据 闽 值 
将 其 分 为 四 类 : 非 适 生生 境 (0 ~ 0.3); REEE 
(0.3 ~0.4); 边缘 适 生 生境 (0.4 ~0.5); WEE 
(0.5 ~0.6) 和 高 适 生生 境 (0.6 ~1)。 






























































































































































































































































2 结果 与 分 析 
2.1 筛选 独立 图 层 并 建立 模型 

TR dis E Bi DSL LER BL AH DR RB Ep BS IR, TE 
95% '& fri EX [8] V3 ffs 3o Eh ii E s [8] 28 vr THE S A JE 
( 表 2) 。 需 要 指出 的 是 ， 有 昼夜 温差 与 年 温差 比值 
代表 了 昼夜 温差 和 年 温差 的 消长 关系 ， 由 于 是 常 
数 ， 可 便于 在 不 同年 际 间 比较 。 在 BIOCLIM, 
DOMAIN 和 MAXENT 物种 分 布 模型 运行 条 件 下 ， 
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对 印度 块 菌 在 云南 省 分 布 现状 进行 模拟 (图 2)， 
并 对 模型 准确 度 进行 对 比 ( 表 3)。 如 果 模 型 曲 
线 下 面积 的 平均 值 越 大 ， 则 说 明 模 型 稳定 性 越 
高 ， 而 模型 的 标准 差 越 小 ， 则 说 明 随 机 取样 对 模 
型 的 影响 程度 越 小 ( 曹 向 锋 等 , 2010) 。 结 果 显 
示 ，MAXENT 模型 在 各 类 模型 组 合 中 都 具有 较 
高 的 拟 合 准确 性 。 对 于 单纯 的 气候 模型 ，MAX- 
ENT 模型 稳定 性 高 于 BIOCLIM 和 DOMAIN， 但 
容易 受到 随机 取样 的 影响 ， 对 于 气候 物种 模型 ， 
MAXENT 模型 获得 了 最 高 的 模型 拟 合 表现 ， 不 
容易 受 随 机 取样 影响 ， 同 时 具有 较 高 的 稳定 性 ， 
并 能 同时 处 理 类 别 变量 和 连续 变量 类 型 ， 故 在 对 


R2 满足 空间 独立 性 的 图 层 


Table 2 The list of spatial-independent layers 








编码 Code 名 称 Name 
AP 年 降水 Annual precipitation 
SLO 坡度 Slope 
SOI 土壤 类 型 Soil types 
PWQ 最 湿 季 度 降 水 Precipitation of wettest quarter 
GEO 地 貌 类 型 Geomorphology 
MTC 最 冷 月 份 最 低温 Min temperature of coldest month 
VEG 植被 类 型 Vegetation types 
MDR 昼夜 温差 月 均值 Mean diurnal range 
ISO 昼夜 温差 与 年 温差 比值 othermality 
UE 最 干 季度 时 平均 温度 


Mean temperature of the driest quarter 





(a) BIOCLIM 





一 一 省 内 各 市 行政 边界 Municipal administrative boundaries 


BS 非 适 生 区 Unsuitable distribution 
EX] 低 适 生 区 Lowly suitable distribution 


(b) DOMAIN 


(c) MAXENT 





边缘 适 生 区 Marginally suitable distribution 








Ex 适 牛 区 Suitable distribution 
EE 高 适 生 区 Highly suitable distribution 





图 2 基准 气候 下 印度 块 茵 

















在 云南 省 潜在 分 布 预测 
Fig.2 Under current climate condition, the predicted potential distribution of T. indium within Yunnan province. A climate model : 


(a) BIOCLIM, (b) DOMAIN, (c) MAXENT; B climate-species model; C climate-topographic model. The habitat suitability 


is presented from red ( highly suitable area) to green (unsuitable area) , standing for various circumstance of habitat 
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印度 块 菌 在 未 来 气候 变化 情景 下 的 空 
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表 3 模型 间 曲 线 下 面积 (AUC) 比较 


Table 3 The evaluation of models accuracy by AUC comparison 














BIOCLIM DOMAIN MAXENT 
模型 组 成 Model components 平均 值 标准 差 平均 值 标准 差 平均 值 标准 差 
Standard Standard Standard 
Mean X Mean E Mean n 
deviations deviations deviations 
气候 模型 Climate model 0. 68 0.029 0. 79 0. 035 0.85 0. 039 
气候 物种 模型 Climate-species model 0. 69 0.029 0.72 0.037 0.95 0.015 
气候 地 形 模型 Climate-topographic model 0.74 0.038 0.65 0. 022 0.93 0.016 


分 布 现状 模拟 及 因子 重要 性 评价 分 析 上 采用 该 模 
型 结果 ; 对 于 气候 地 形 模 型 ，MAXENT 模型 表 
现 略 差 于 气候 物种 模型 模拟 ， 但 模拟 结果 也 不 易 
受 随 机 取样 影响 。 根 据 拟 合 程度 最 优 的 模型 来 
看 ， 预 测 结果 显示 云南 省 中 部 地 区 是 主要 的 适 生 
区 域 , 集中 于 迪 庆 藏族 自治 州 ， 怒江 傈 傈 族 自 治 
州 ， 大 理 白 族 自治 州 ， 楚 雄 莱 族 自治 州 ， 丽 江 市 ， 
昭通 市 ， 曲 靖 市 ， 昆 明 市 ， 玉 溪 市 和 保山 市 。 其 
中 ， 高 度 适 生 区 体现 了 最 佳 的 环境 因子 及 最 适宜 
的 气候 因子 组 合 ， 在 云南 省 呈 高 度 分 散 分 布 。 
2.2 ”因子 重要 性 评价 

MAXENT 模型 自 带 的 Jackknife 功能 可 辅助 
分 析 因 子 重要 性 〈 图 3) ， 结 果 显 示 了 当 剔 除 因 
子 、 保 留 因子 及 保留 全 部 因子 时 对 模型 的 增益 值 
(Gain) 及 模型 准确 度 (AUC) 的 贡献 ,综合 
虑 以 上 两 个 方面 ， 年 降水 、 最 湿 季 度 降 水 、 最 冷 
月 份 最 低温 、 地 和 瑶 类 型 及 土壤 类 型 具有 最 大 影响 
力 ， 而 昼夜 温差 与 年 温差 比值 则 作用 最 小 。 和 其 
它 因 子 不 同 ， 当 剔除 坡度 因子 时 ， 模 型 总 的 准确 












































度 上 升 ， 故 坡度 因子 的 加 入 虽然 能 增加 模型 增益 
值 ， 但 同时 也 减少 了 最 终 准确 性 。 然 而 ， 根 据 专 
家 知识 和 野外 采集 经 验 ， 对 印度 块 菌 影响 较 大 的 
因子 除了 气候 、 地 形 地 貌 因子 外 还 应 包含 植被 因 
子 ， 而 在 模型 结果 中 未 能 体现 ， 其 原因 可 能 归 为 
以 下 两 点 : (1) 模型 输入 的 云南 省 植被 图 制作 
时 间 久 远 ， 难 以 将 目前 印度 块 菌 共生 植被 的 情况 
正确 展现 ; (2) 在 大 尺度 上 的 空间 模拟 ， 植 被 
同样 受 控 于 气候 和 地 形 因子 限制 ， 故 可 能 导致 相 
对 重要 性 的 减 小 。 同 时 ， 除 了 分 析 各 因子 对 模型 
总 体 增益 值 和 准确 度 的 影响 ，MAXENT 模型 还 

能 够 进一步 分 析 类 别 变 量 (如 植被 类 型 、 土 壤 
类 型 等 ) 对 印度 块 菌 发 生 概 率 的 影响 ， 如 果 最 
终 概 率 达到 50% 以 上 ， 则 说 明 该 类 别 是 印度 块 
菌 生 长 和 分 布 的 适宜 类 型 。 模 型 模拟 结果 显示 针 
叶 林 和 常 绿 阔 叶 林 较 其 它 植被 类 型 更 适合 印度 块 
菌 生 长 ， 同 时 中 山 和 喀斯特 地 形 为 适宜 地 貌 类 
型 ， 紫 色 土 和 棕 壤 为 适宜 土壤 类 型 。 

2.3 未 来 气候 情景 下 的 印度 块 菌 空 间 分 布 模式 
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Regularized training gain 





8 剔除 变量 Without variable 


E 保留 变量 With only variable 


0.45 0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85 
曲线 下 面积 
AUC 


8 所 有 变量 With all variables 


图 3 基准 气候 下 ，MAXENT 模型 中 各 指示 因素 评价 


Fig.3 Evaluation of environmental factors in MAXENT model 


a) factors’ contribution to regularized training gain; b) factors" influence on Area Under Curve ( AUC) 
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在 未 来 气候 变化 情景 下 ， 通 过 模拟 从 基准 Mn 境 类 型 ，A2 情景 呈 组 慢 递 减 趋势 ， 而 


线 气 候 到 2080 年 印度 块 菌 分 布 模式 (图 4)，A2 
和 B2 情景 显示 了 未 来 印度 块 菌 新 增生 境 主要 呈 





2 情景 则 先 缓慢 减少 ， 
B2 ' 





北上 趋势 ， 更 加 的 边缘 化 和 分 散 化 ， 临 沧 市 、 普 
洱 市 的 局 部 地 区 可 能 成 为 未 来 的 潜在 生境 区 。 在 
A2 气候 情景 下 ， 减 少 的 生境 将 主要 出 现在 丽江 


市 、 昆 明 市 和 玉溪 市 ;而 在 B2 气候 情景 下 
少 的 生境 集中 于 丽江 市 、 大 理 白 族 自治 州 、 


玉溪 


市 和 曲靖 市 。 为 了 进一步 了 解 生境 面积 变化 的 具 
体 情况 ， 在 A2 和 B2 气候 情景 下 ， 计 算出 每 个 时 
段 生 境 占 总 区 域 的 面积 百分比 〈 表 4) ， 对 于 未 





基准 气候 -2080 年 
Baseline-2080 





A2 


图 4 ， 从 基准 气候 条 件 至 2080 年 印度 块 菌 9 
Fig.4 The habitat shifting pattern of future T. indium distribution. From baseline to 2080 


后 快速 增加 ; 对 于 非 适 
青 景 减少 的 速度 高 于 A2; 对 于 未 来 无 


变化 区 域 ，B2 o A2; 对 于 未 
青 景 较 为 稳定 ， 没 有 太 大 


来 新 增生 境 类 型 ， 
变化 ， 而 B2 情景 


apap. 呈 先 增加 后 减少 


趋势 。 而 通过 分 布 概率 建立 的 印度 块 菌 未 来 分 布 
分 比 数据 ( 表 5)， 
在 A2 情景 下 ， 虽 然 每 个 时 段 增长 速率 有 波动 ， 


预测 图 (图 5) 及 相应 的 面积 


生境 类 


型 总 


ED 


体 偏 向 于 边缘 适 生 区 和 适 生 区 类 型 ， 





而 在 B2 情景 下 ， 生 境 适 生 程度 总 体 降低 。 


B2 








省 内 各 市 行政 边界 


Municipal administrative boundaries 


未 来 减少 生境 


Future decreased distribution 








非 适 生 区 


Unsuitable distribution 


pa 未 来 无 变化 生境 


Future inalterable distribution 


ma 术 米 新 增生 境 


Future increased distribution 


E 境 变化 模式 预测 图 ( 从 二 值 图 层 推 导 ) 


under A2 and B2 scenario (Calculated from binary maps) 


表 4 各 生境 类 型 面积 所 占 总 研究 区 域 百分比 


Table 4 Different habitat percentage of whole study area 


生境 类 型 Habitat types 


A2 


B2 


























Baseline-2020 2020-2050 2050-2080 Baseline-2020 2020-2050 2050-2080 
未 来 减少 生境 Future decreased distribution 3.28 2.11 1.01 1.62 0. 90 3.40 
非 适 生 区 Unsuitable distribution 82.39 82.67 81.65 82.11 79.44 78.35 
未 来 无 变化 生境 Future unaltered distribution 11.12 12.21 14.21 12.79 15.36 16.26 
未 来 新 增生 境 Future increased distribution 3.20 3.00 3.14 3.48 4.30 2.00 
表 5 未 来 气候 情景 下 适 生生 境 面积 百分比 变化 
Table 5 The suitable habitat area percentage changes under different climate change scenarios 

适 生 区 类 型 Habitat suitability Baseline p 

2020 2050 2080 2020 2050 2080 
非 适 生 区 Unsuitable distribution 74.01 78.73 74.61 74.41 73. 85 74. 17 76.31 
低 适 生 区 Lowly suitable distribution 9.07 6.87 8.55 8.42 8.32 8.54 7.57 
边缘 适 生 区 Marginal suitable distribution 6.38 5.42 6.49 6.73 6. 66 6.62 5.99 
3 ^E [X Suitable distribution 4.45 3.90 4.81 4.74 5.15 4.94 4.59 
高 适 生 区 Highly suitable distribution 6.09 5.08 5.55 5.70 6.02 5.74 5.48 
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图 5 在 未 来 不 同 气候 情景 下 印度 块 菌 在 云南 省 潜在 分 布 预测 图 (从 分 布 概率 推导 ) 


Fig.5 Under future climate change scenarios, the predicted potential distribution of T. indium within Yunnan province (Based on 


distribution probability). A Scenario A2; Baseline-2020, 2020-2050, 2050-2080; B Scenario B2. The habitat suitability 


is presented from red ( highly suitable area) to green (unsuitable area) , standing for various circumstance of habitat 


3 讨论 

基于 本 研究 建 模 前 所 提出 的 3 点 假设 ， 发 现 
印度 块 菌 未 来 新 增生 境 具 有 高 海拔 分 布 趋势 ， 有 
部 分 研究 将 其 它 物种 的 类 似 现象 归 为 对 气候 变化 的 
响应 机 制 (Mideley 等 ，2002; Thomas ^&, 2004) 。 
不 同时 段 的 印度 块 菌 ， 其 适 生 区 和 高 适 生 区 域 与 
基准 线 气 候 条 件 下 的 分 布 区 范围 差异 较 小 ， 是 未 
来 无 变化 生境 的 主要 组 成 部 分 。 体 现 了 印度 块 菌 
对 环境 变化 的 耐 受 能 力 较 强 ， 呈 现 积 极 适应 并 扩 
张 原 有 自然 生境 的 趋势 。 同 时 ， 需 要 明确 的 是 ， 
正 是 由 于 建 模 前 所 具有 的 假设 条 件 ， 故 模型 体现 
的 只 是 未 来 印度 块 菌 可 能 分 布 的 模式 之 一 。 在 物 
种 实际 扩散 途中 ， 可 能 会 因为 受到 人 为 干扰 、 自 
然 灾 害 、 病 虫害 、 竞 争 物种 等 干扰 ， 而 阻碍 其 在 
适宜 生境 区 的 生长 ， 因 此 ， 采 纳 模型 预测 结果 时 
须 持 慎重 态度 。 同 时 ， 未 来 印度 块 菌 生境 发 生 向 
北 迁 移 ， 可 能 会 对 该 区 域 原 有 植 生 共生 菌 也 产生 
影响 。 如 Suzuki (2002) 的 研究 表明 真菌 群落 和 

















植物 群落 类 似 ， 具 有 宏观 和 微观 太 度 上 的 演 替 。 
外 生 菌 根 真 菌 需要 从 共生 树种 获取 营养 来 源 ， 故 
共生 树种 的 演 蔡 对 真菌 的 演 蔡 起 到 重要 作用 ， 也 
对 真菌 间 的 竞争 产生 促进 或 抑制 作用 (Prentice, 
1992)。 因 为 具体 讨论 未 来 可 能 的 真菌 演 蔡 现象 
已 超出 了 本 文 所 能 回答 的 范围 ， 故 需要 借助 于 今 
后 的 相关 研究 才能 解答 。 

虽然 国内 外 对 主要 动 植物 物种 的 生境 建 模 有 
很 多 尝试 (Anderson 等 , 2002; Thuiller 等 ,2003 ; 
Leathwick 等 ，2005; Pearson 等 ，2007; Cayuela 
等 , 2009) ， 但 很 少 有 研究 运用 地 理 信 息 系 统 和 
物种 分 布 模型 来 模拟 林 下 资源 的 潜在 生境 ， 原 因 
来 自 于 多 个 方面 ， 如 基础 数据 不 全 、 对 物种 认识 
不 全 及 生境 的 特殊 性 等 ( Yang 等 , 2006)。 对 于 
林 下 植物 类 型 ， 难 以 直接 从 遥感 影像 上 实现 不 同 
时 段 和 空间 的 监测 ， 而 传统 的 生态 调查 数据 通常 
建立 在 不 同 研究 区 域 ， 在 不 同时 期 内 由 不 同 采集 
者 收集 ， 在 时 空 上 不 具有 连续 性 。 在 更 多 情况 
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下 ， 趋 于 采集 的 困难 性 ， 研 究 者 对 于 感 兴趣 的 研 
究 对 象 难以 得 到 足够 空间 分 布 的 信息 ， 从 而 导致 
了 资源 的 利用 和 管 护 不 当 ， 这 是 国内 外 所 面临 的 
共同 问题 。 

在 建立 模型 过 程 中 ， 需 要 就 各 种 模型 的 优点 
和 局 限 性 进行 优化 ， 必 要 时 考虑 模型 组 合 ， 以 期 
提高 预测 能 力 。 本 研究 最 终 第 选 的 MAXENT fi 
型 ， 能 处 理 较 小 样本 集 ， 在 同时 分 析 类 别 变量 和 
连续 变量 时 能 获取 较 高 拟 合 效果 ， 但 要 求 独 立 样 
本 数据 以 及 独立 的 气候 变量 。 不 平衡 的 样本 效应 
以 及 高 度 自 相关 的 气候 变量 将 导致 无 效 模型 
( Jiménez-Valverde 和 Lobo, 2006; Loiselle 等 ， 
2008), ， 这 一 问题 可 以 通过 改进 训练 数据 集 的 质 
量 和 提取 独立 性 的 环境 变量 来 完成 ， 在 这 种 情况 
下 ， 对 实验 对 象 分 布 的 认 知 程度 和 专家 知识 非常 
重要 。 同 时 ， 如 果 输 入 数据 图 层 的 制作 单位 、 制 
作 时 间 、 空 间 分 辨 率 差 异 较 大 ， 也 会 给 模型 引入 
很 大 的 不 确定 性 。 将 模型 模拟 结果 与 野外 实测 结 
果 进 行 比 对 是 对 模型 准确 性 验证 的 方式 之 一 ， 实 
验 室 建 模 与 野外 实测 过 程 是 密切 联系 ， 环 环 相 扣 
的 。 此 外 ， 通 过 建立 地 理 信息 数据 库 ， 可 以 实现 
不 同 区 域 、 地 方 的 数据 共享 ， 能 更 为 全 面 的 了 解 
物种 分 布 与 扩散 信息 。 而 模拟 未 来 情景 下 的 物种 
分 布 模式 ， 对 今后 开展 保护 工作 及 制定 相关 决策 
具有 重要 意义 。 























4 结论 

(1) 在 对 印度 块 菌 生境 分 布 现 状 模拟 过 程 
中 ，MAXENT 模型 具有 最 优 的 表现 ， 模 型 稳定 
性 高 且 不 易 受 随机 取样 影响 。 当 添加 了 其 它 环境 
因子 后 ， 模 型 的 总 体 预测 能 力 提高 。 在 今后 的 研 
究 中 ， 除 了 和 气候 条 件 ， 物 种 间 的 相互 作用 也 是 建 
模 过 程 中 需要 考虑 的 方面 。 

(2) 在 对 模型 模拟 结果 进行 分 析 时 ， 需 要 
同时 依靠 专家 知识 和 野外 工作 经 验 ， 并 就 不 断 变 
化 的 环境 因子 ， 如 植被 类 型 ， 不 断 更 新 以 提高 模 
型 预测 力 。 本 研究 模拟 结果 显示 在 大 的 空间 尺度 
上 ， 年 降水 、 最 湿 季度 降水 、 最 冷 月 份 最 低温 、 
地 貌 类 型 及 土壤 类 型 具有 最 大 影响 力 。 

(3) 在 A2 和 B2 情景 下 ， 印 度 块 菌 未 来 分 
布 将 呈 北 上 趋势 ， 同 时 新 增生 境 更 为 边缘 化 和 分 
散 化 。 两 类 气候 情景 下 ， 生 境 类 型 的 变化 速率 不 




















同 ， 总 体 上 看 ，A2 气候 情景 较为 稳定 ， 未 来 生 
境 类 型 中 边缘 适 生 区 和 适 生 区 类 型 比例 增加 ， 而 
B2 气候 情景 下 生境 面积 变化 相对 较 大 ， 同 时 生 
境 的 适 生 程度 减低 。 
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